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u _ 

study the effect of the matter term on the evolution of the solar neutrinos 
when the neutrino parameters are those of the 'just-so' case. The extreme non- 

(N ■ 

adiabatic effects at the edge of the sun reduce the expression for the survival 

T 1 

' probability in the just-so case to that of the vacuum case. This conclusion 

o ' 

CN ] is independent of the width of the extreme non-adiabatic region, which is a 

' 

o 
o 
o 



function of the density profile of the sun beyond r > 0.9Rs- However, in 
its propagation, neutrino encounters regions of moderate (non-extreme) non- 
adiabticity. Neutrino traversal through these regions give corrections to the 
survival probability which are profile dependent. 



I. INTRODUCTION 

The solar neutrino data of Super-Kamiokande has lead to a reexamination of all the 
solutions of the solar neutrino problem In addition to the overall suppression rate, Super- 
Kamiokande has measured the day-night asymmetry and the spectrum of the scattered 
electrons. These two measurements are independent of the overall normalization of the 
neutrino flux. This normalization has the largest uncertainty among all the predictions of 
the solar models 0J^. Hence data, which are independent of it, are extremely important in 
studying the properties of neutrinos from the sun. The present data, especially the recoil 
spectrum, favor 'just-so' oscillation solution to the solar neutrino problem 



In analyzing the solar neutrino data in terms of 'just-so' oscillations, the expression 
used for electron neutrino survival probability Pee is simply the vacuum survival probability. 
Given that the matter term in the sun is several orders of magnitude greater than the 
mass-squared difference, the question arises whether the use of vacuum survival probability 
is justified. A very simple justification can be given the following way: The mass-squared 
difference 6 required for just-so oscillations is about 10^^° eV^. The matter term 



is much larger than S as long as the density p is greater than 0.01 gm/cc. First we assume 
that the density is spherically symmetric and it falls abruptly to once it decreases to a 
value of 0.01 gm/cc. As long as the neutrino is in the sun, A completely overwhelms 6 and 
the electron neutrino in the sun is essentially the higher mass eigenstate. Hence during its 
travel through the sun, the neutrino does not oscillate and it simply acquires a phase. Thus 
throughout the travel through the sun, an electron neutrino remains an electron neutrino. 
When the neutrino comes out of the sun, it passes through the abrupt change in density 
from 0.01 to 0. In such a case, the flavour composition is unchanged. Since the neutrino in 
the sun remained an electron neutrino, we have an electron neutrino coming out of the sun. 
Thus the starting point of neutrino evolution is transferred from the core of the sun to the 
edge of the sun 0. 

But now we can raise the question: Suppose the solar matter density does not change 
abruptly to zero but goes to zero smoothly. Then at some point in its propagation the 
neutrino will pass through a region where the density is equal to the resonant density. 
In such a situation, how is the neutrino oscillation probability modified? We address the 
question below. 



For simplicity here we consider only two flavor oscillations. The flavor states z/g and z/^ 
are linear combinations of the two mass eigenstates ui and z/2 
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(1) 



II. JUST-SO SURVIVAL PROBABILITY 



— COS 9 sin 9 V2 

i/^ = — sin 9 vi + cos 9 U2. 



(2) 



Without loss of generality we can take V2 to be more massive than Vi and hence 5 = 
ml — ml > 0. If ^ < 7r/4, then Ue is predominantly the lighter state. If ^ > 7r/4, then i/e 
is predominantly the heavier state. Thus the two physically distinguishable possibilities are 
covered by taking the range of 9 to be (0, 7i/2). Since we are considering 'just-so' oscillations, 
we will assume that S = 10~^° eV^. We will also assume that the value of 9 is moderately 
large, i.e. 9 is not finetuned to be very small or be very close to 7r/4. 

Electron neutrinos are produced in the core of the sun. Since the solar matter near the 
core is very dense, one must take the matter term A into account in determining the mass 
eigenstates in the core. We can define instantaneous matter dependent mass eigenstates 



Ue = COS 9m Vim + Siu 9m l^2m 
1^1^ = - sin 9m l^lm + COS 9m l^2m, 

where 9m is the matter dependent mixing angle and is given by 

5cos29-A 
cos 29m = r ■ 

Om 

Sm is the matter dependent mass-square difference and is given by 

Sm = J{Scos29 -Ay + {6sm29y. 
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As the neutrino travels through the sun, it encounters matter of constantly decreasing den- 
sity. In the basis of matter dependent mass eigenstates, the evolution equation for the 
neutrino is 
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-6m{t) -AiE9m{t) Vim 
UE9m{t) 5m{t) V2m 

The propagation of neutrino is adiabatic as long as the off-diagonal terms are smaller than 
the difference of the diagonal terms in the above equation, i. e. 



2EA5sin2e 

Om > - 



52 

In the body of the sun the density profile is of the form 



(7) 



p(r) =poexp(-10.54r/i?,) (8) 

where po is the core density (about 140 gm/cc) and Rg is the radius of the sun (about 7 x 10^ 
km) Therefore 

10.54 
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'a 




P 



3 X 10"'^ eV. (9) 



As long as p > 0.01 gm/cc (or r < 0.9i?s) the matter term A is much larger than 6, hence 
6m — A. According to our assumption on 9, sm29 ~ 1. Substituting all these in Eq. (^, 
the adiabatic condition becomes 

> 6 X 10"^^(E/MeV). (10) 

The most compelling reason for considering just-so oscillations as a solution to the solar 
neutrino problem is Super-Kamiokande data on electron spectrum due to neutrinos. 
The observed range of these neutrinos is 6-14 MeV and they are peaked around 10 MeV. 
Henceforth, we will take E = 10 MeV for illustrative purpose. Conservatively we require 
that A"^ should be at least 10 times the value on RHS. With this requirement, we find that 
the adiabatic condition is satisfied if 

p > 0.01 gm/cc. (11) 

The adiabatic condition is likely to break down if p < 0.01 gm/cc, or for radial distances 
greater than r > O.dRg. 

When the neutrino is well out of the sun the matter term is much smaller than 6 and 
adiabaticity is restored. In such a case we have, essentially, vacuum propagation. The radial 
distance at which this happens is a function of the density profile of the sun beyond O.dRs- 
For example, if the density falls linearly beyond 0.9i?s, vacuum propagation starts around 



r = Rg. If the density fall continues to be the exponential form shown in Eq. then 
vacuum propagation starts only around 2Rs. Between the time of breakdown of the adia- 
baticity and the time of its restoration, the off-diagonal terms in the equation of propagation, 
Eq. (^), are comparable to the diagonal terms. In fact, for some intermediate region, the 
off-diagonal terms completely dominate, making the propagation extremely non-adiabatic. 
These extreme non-adiabatic effects cause the just-so oscillation probability to reduce to 
vacuum oscillation probability. 

An electron neutrino is produced in the core of the sun at time to and proipagates 
adiabatically upto ti. Between times ti and t2 its propagation is non-adiabatic. Beyond t2 
the neutrino propagates in vacuum and is detected at ^3. Its state vector at ^3 is 

I*e(t3)) = E \^j)('^P H^Ah - t2)) M,,exp (-1 r 6f{t)dt) Ug, (12) 

hi ^ ^ 

where sf{t) are the matter dependent energy eigenvalues in the sun and Si are the energy 
eigenvalues in vacuum. Mjj is the amplitude for non-adiabatic transition from mass eigen- 
state z/jm at t = ti to mass eigenstate Vj at t = t2- Ug is the e — i element of the matter 
dependent mixing matrix at the core of the sun. From the above equation, we obtain the 
electron neutrino survival probability to be 

Pee=\{l^e\^e{h))\'' 

= E E UejMijUgUefM,,j,Ugexp |^ {e^t) - ef{t))dt^ | exp {^(e,. - e,) (ta - ts)} • (13) 

ji 

Because of the variation in the production region, the phase picked up in the time from to 
to ti can be averaged out. So exp (i Jt^{£it{t) — ef{t))dtj can be replaced by 6^. Thus, Pee 
is simplified to 

Pee = E UejM,,UgUefM,^,Ugexp {l{Ej, - ej){t^ - h)} . (14) 

The mixing angle at the core of the sun 6c = 7r/2, because the electron neutrino is the 
heavier mass eigenstate due to the dominance of the matter term. In vacuum, of course, it 
is 9. Substituting these in the expression for Pee, we get 
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Pee = Xi2 cos^ 9 + {1 — X12) sin^ 9 

+ 2 sin 9 cos 9^ (Mi*2M22) COS (^2.54^ j , (15) 

where xu = |Mi2p is the probabihty for the I'lmiti) evolve into z/2 at t2- In the above 
equation we have also made the extreme non-relativistic apporoximation for neutrinos and 
replaced the difference in energies by the mass-squared difference and the time time of travel 
(^3 — h) by the distance of travel x. 

The survival probability given in Eq. (|15|) should give us vacuum survival probability 
if the density falls abruptly to 0, that is t2 ~ is extremely small. In such a case, it was 
shown that M12 = {UlU)i2 = sin(^^5 — 9) = cos9 because 9s, the mixing angle in the sun, 
is 7r/2 [0]. For M12 = cos^^. Pee immediately reduces to vacuum survival probability. Now 
the question is: what is M12 if the density falls smoothly to zero and the distance between 
the point of breakdown of adiabaticity and the point of start of vacuum propagation is a 
significant fraction of solar radius. In such a case, the calculation of M12 is complicated and 
depends on the density profile of the sun beyond 0.9i?s. Let us calculate M12 in a simplified 
situation where we will assume that the off-diagonal terms are much larger than the diagonal 
terms between the point of breakdown of adiabaticity (ti) and the point of restoration of 
adiabaticity {12). That is, in the region where adiabatic approximation is not valid, the 
evolution is extremely non-adiabatic. Then the evolution equation is 
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(16) 



Integrating this equation from ti to ^2, we get 

Mi2 = sin(^(ti)-^(t2)). (17) 

As mentioned above, until the time ti, A » 5, hence 9{ti) = 7r/2. At time t2, vacuum 
propagation starts, hence 9{t2) = 9. So we get M12 = cos6', which reduces Pee to vacuum 
survival probability. Note that we have not made any assumption about how the density 
varies to zero. The width of non-adiabatic region (given by (ti — ^2)) is irrelevant. In 
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case 9 < Ti/A there may be a point within this non-adiabatic region where the resonance 
condition is satisfied. However, the discussion above is completely independent of whether a 
resonance exists or not. The reason for this is that the breakdown of adiabaticity condition 
is independent of the existence of resonance. This is in contrast to the usual MSW effect 
where the adiabaticity condition breaks down only in the neighbourhood of a resonance. 

The above simplication is not very good because, just after the breakdown of adiabatic- 
ity and also just before the restoration of the adiabaticity, the diagonal and off-diagonal 
terms in the evolution matrix are comparable. We need to consider these regions of 'mod- 
erate non-adiabaticty' separately. To study how good the above assumption of extreme 
non-adiabaticity throughout the region of non-adiabaticity is, we define a quantity Non- 
Adiabaticity Paramater (NAP). It is the ratio of the non-diagonal term to the diagonal 
term in the evolution matrix (Eq. 



NAP 

-■ (18) 
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m 

Note that NAP is proportional to dp/dt and is much smaller than 1 for regions of adia- 
baticity and is comparable to or greater than 1 in the region where adiabatic approximation 
breaks down. The region where NAP is much greater than 1 is the region of extereme 
non-adiabaticity. Suppose density profile varies smoothly upto the edge of the sun and then 
abruptly falls to zero. If the adiabatic approximation holds till the edge of the sun, the 
graph of NAP vs r looks like a Dirac delta function aX r = Rg. In the discussion below, we 
consider two different density profiles. 

• Exponential density fall (given in Eq. (|D) for all r, even beyond r = 0.97?^. 

• Exponential density fall upto r = 0.9-Rs and then a linear density fall beyond r = 0.97?^. 
The equation of the linear fall is determined by requiring that the value of the density 
and its first derivative should match at r = 0.9/?,. 
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We see that for linear fall, NAP is zero until r c:^ Rg and then rises sharply at r = Rs, which 
resembles the NAP graph for an abrupt density change. Hence for linear fall at the edge 
of the sun, the just-so oscillation probability reduces to the vacuum oscillation probability. 
The region of moderate non-adiabaticity is very narrow (r = 0.957 Rg to r = O.99I-R5). 

However for exponential fall, NAP rises slowly, becomes greater than 0.3 (breakdown of 
adiabatic approximation) at r = 0.96Rs, crosses 10 around r = 1.12Rs and falls below 10 
again for r = 1.71-Rs- The adiabatic approximation (NAP j 0.3) is restored only for r = 
2.04i?s. We see that regions of moderate non-adiabaticity are reasonably large both during 
increasing NAP (O.lGi?^) and during decreasing NAP 0.33-Rs- The transition amplitude 
between the two edges of extreme non-adiabaticity can be easily calculated, as shown in 
Eqs. ([T6| ) and (jl^)- However, the effects of propagation through the regions of moderate 
non-adiabaticity (0.3 < NAP < 10) must be taken into account. Since the regions of 
moderate non-adiabaticity are wide, these can be calculataed only by numerical integrattion 
of the evolution equation in these regions. These moderate non-adiabatic effects can lead 
to the energy dependence of the just-so survival probability to be different from that of the 
vacuum survival probability. This investigation is in progress. 

As mentioned above, for exponential density fall for all r, the survival probability for 
just-so oscillations is likely to be different from that of vacuum oscillations. Recently a 
calculation of such a modified survival probability was presented P] . We find the calculation 
in P] untenable for the following reasons. 

1. The calculation uses an expression for non-adiabatic jump probability which was cal- 
culated under the assumption that the physical region in the sun of extreme non- 
adiabaticity is very small so that it can be taken to be a point (namely the resonance 
point) and the non-adiabatic jump is calculated only at this point. However, as can 
be seen from Figure 1, the physical region of non-adiabaticity is of the order of solar 
radius, hence the assumption of non-adiabatic jump occuring at a point is unreason- 
able. 
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2. In ref it was pointed out that the survival probabihty is asymmetric under the 
exchange 9 ^ -n jl — Q. The expression for survival probability uses a standard formula 
for non-adiabatic jump probability. However, this formula is derived assuming the 
existance of a resonance Hence it cannot be used for scenarios where Q > it/ 4. 

III. CONCLUSIONS 

• If the solar density abruptly falls to zero at the edge of the sun, just-so oscillation 
probability reduces to vacuum oscillation probability. 

• If the density fall is smooth, the just-so oscillation probability reduces to vacuum 
oscillation probability if one assumes that the evolution is extremely non-adiabatic 
throughout the region where adiabatic condition is not valid. This conclusion is inde- 
pendent of the exact density profile near the edge of the sun or the exact width of the 
region of non-adiabaticity. 

• The regions of moderate non-adiabaticity can give corrections to the simple vacuum 
oscillation probability. However, for linear density fall these corrections are negligible. 

• If the recoil spectrum analysis of Super Kamiokande points towards a departure from 
the simple vacuum oscillation probability for just-so parameters, these departures can 
be used to obtain information on the density profile at the edge of the sun. 

Acknowledgements: We thank Prof. Paul Langacker for a helpful discussion. 
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FIG. 1. Plot of NAP X {= r/Rs), for linear density fall (solid line) and exponential density 
fall (dotted line). 
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